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Abstract 



In this paper we present the master integrals necessary for the analytic 
calculation of the box diagrams with one electron loop (Np = 1) entering in 
the 2-loop (a 3 ) QED virtual corrections to the Bhabha scattering amplitude 
of the electron. We consider on-shell electrons and positrons of finite mass 
■ m, arbitrary squared cm. energy s, and momentum transfer t; both UV 

and soft IR divergences are regulated within the continuous D-dimensional 
regularization scheme. After a brief overview of the method employed in the 
calculation, we give the results, for s and t in the Euclidean region, in terms 
of 1- and 2-dimensional harmonic poly logarithms, of maximum weight 3. The 
corresponding results in the physical region can be recovered by analytical 
continuation. For completeness, we also provide the analytic expression of 
the 1-loop scalar box diagram including the first order in {D — 4) . 
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1 Introduction 



The Bhabha scattering process plays a key role in the phenomenology of particle 
physics, since it is employed in order to determine the luminosity of electron-positron 
colliders. In particular, the small angle Bhabha scattering has been used to measure 
the luminosity at high energy colliders, such as LEP and SLC The large angle 
Bhabha scattering is instead employed in measuring the luminosity of flavor factories 
(BABAR, BELLE, BEPC/BES, DA$NE, VEPP-2M). 

The accuracy of the luminosity measurements depends on the precision of the 
theoretical predictions for the Bhabha scattering cross section, since the luminosity 
is defined as the ratio of the number of events observed and the theoretical cross 
section for the Bhabha process. 

For this reason, in the last three decades, a number of publications have been 
devoted to the study of the radiative corrections to this process. At the one-loop 
level, the complete set of radiative corrections to the Bhabha scattering cross sec- 
tion, in the framework of the Standard Model of the electroweak interactions, has 
been calculated several years ago jT]; at the 2- loop level, some results, namely the 
factorisable subset [2], have been obtained, but a complete and exact evaluation of 
the 2-loop QED quantum corrections is still missing. 

Within the approximation of only keeping contributions enhanced by factors of 
\rt(s/m 2 e ) a large amount of work has already been done. In this approximation 
the 2-loop corrections to the large angle Bhabha scattering cross-section containing 
electron-positron pairs were considered in [Hill]. These calculations include correc- 
tions coming from the interference of the graph in Fig. ^with the tree level diagrams. 
The real hard-pair production was studied in [5|. The contributions without pair 
production, virtual or real, were considered in Ej, where, however, the 2-loop 
double box graphs (planar and crossed) were ignored. 

Recently, the full set of 2-loop virtual QED corrections to Bhabha scattering 
has been calculated, but with the approximation of neglecting the electron mass jH] . 
The second order logarithmic corrections to the large angle Bhabha scattering cross 
section, coming from graphs that do not involve vacuum polarization insertions, 
have been calculated in j^j. 

The 2-loop 4-point Feynman diagrams are one of the essential ingredients in 
the calculation of the next-to-next-to-leading order QED corrections to the Bhabha 
cross section; not surprisingly, their evaluation represents one of the main tech- 
nical difficulties encountered in calculating these corrections. With the exception 
of cases presenting specific kinematic configurations, the calculation of 2-loop box 
diagrams is an open problem. In the recent past, the complete evaluation of the 
master integrals for the 2-loop 4-point functions, with massless propagators, has 
been carried out; this has been done, analytically, in the case of massless external 
legs [TOJ HH C21 EE] , and in the case of three massless and one off-shell external legs 
^HlCSl- Numerical results, in the non-physical region s, t < 0, were also presented 
in [T3|. For what concerns 4-point diagrams with massive propagators, as far as 
the authors know, the only available result is in [17] . where the scalar double box 
integral with four massive and three massless internal lines, and external legs on 
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Figure 1: 2-loop box diagram relevant for the purpose of the paper; the momenta 
Pi and p 2 are incoming, p 3 and p 4 outgoing. 



their mass-shell, is evaluated in terms of Nielsen's polylogarithms (of non-simple 
arguments and maximum weight 3). 

In the present paper, we evaluate the master integrals (Mis) necessary for the 
calculation of the box diagrams with one closed electron loop (Np = 1) entering the 
2-loop virtual corrections to the electron Bhabha scattering amplitude in QED; the 
calculation is carried out without neglecting the electron mass m. Performing the 
calculation without considering the electron massless allows to control the collinear 
singularities. 

The relevant t-channel Feynman diagram is shown on Fig. ^ (the remaining t- 
channel diagrams and the s-channel diagrams can be recovered by crossing), where 
we consider the scattering of an incoming electron of momentum p\ and a positron 
of momentum p 2 , into an outgoing electron and a positron of momenta p% and p4, 
respectively. All the external legs are on their mass-shell, p 2 = —m 2 ; we further 
define 

P = Pi+p 2 , Q=Pi-Ps, s = -P 2 , t = -Q 2 . (1) 

We carry out our calculation in the non-physical region P 2 ,Q 2 > 0; the physical 
region for the Bhabha scattering, s > 4m 2 , t < is to be recovered by analytical 
continuation. 

The interference of this class of diagrams with the tree-level amplitude (which 
provides the 0(a 4 ) contribution to the cross-section we are interested in) can be 
expressed in terms of a large number of 2-loop scalar integrals associated with the 
considered graphs. Following a by now standard approach, we express all the scalar 
integrals that appear in the problem as combinations of a small number of inde- 
pendent scalar integrals, the so-called Master Integrals (Mis) of the diagrams under 
consideration. The reduction procedure that allows to express the generic scalar 
integral in terms of Mis has been discussed extensively in [TT] ITS] , and it is based 
on the use of the Integration by Parts Identities (IBPs) [15] . the Lorentz Invari- 
ance Identities (LI) [TI], and the symmetry properties [TB] of the scalar integrals 
encountered in the problem. The analytic calculation of the Mis is then performed 
by means of the Differential Equations Method [2U1 I2TI 1221 ITT] . 

All the integrals considered in this work are Euclidean, regularized within the 
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dimensional regularization scheme j2H], in which both UV and IR divergences are 
regulated by the same parameter D, the (continuous) number of space-time dimen- 
sions. The results are given as a Laurent series in (D — 4), and the coefficients of 
these series are expressed in terms of generalized 1- and 2-dimensional harmonic 
polylogarithms (HPLs) j2H 1213 12H1 H3 12H| , a suitable generalization of the Nielsen's 
polylogarithms [2H1 EH EH E2]- As we work in the unphysical Euclidean region, 
where s < 0, all the integrals are real. 

The present paper is structured as follows: in Section El we review briefly the 
procedure that allows to express a scalar integral in terms of the Mis, focusing 
our attention on the specific case under consideration; in Section El we review the 
method of differential equations for the calculation of the Mis, giving an explicit 
example in Subsection 13. 11 where the solution of the system of differential equations 
for a 5-denominator four-point function is discussed. In Section E] we provide the 
expression of the 6-denominator scalar integral associated to the (unrenormalized) 
graph of Fig. ^ Appendix [X] contains the definitions of the propagators in terms of 
the loop momenta; in Appendix [B] we give the expression of the 1-loop QED box 
scalar diagram (two massless and two massive internal lines) up to the first order in 
(D — 4) included. Finally, in Appendix we briefly review the formalism of 1- and 
2-dimensional harmonic polylogarithms. 

The complete expression of the contributions of the two-loop box diagrams to 
the corrections of 0(a 3 ) to the Bhabha scattering amplitude will be given elsewhere. 

2 The reduction to the Master Integrals 

In this section, we give a brief overview of the reduction procedure that has been 
employed in order to express, in terms of Mis, the scalar integrals involved in the 
calculation of the a 3 contributions to the Bhabha scattering amplitude. This topic 
is discussed in greater detail in jlT] IIB]. 

At first, we remind the reader that, according to the definition given in [IB] , 
there is just one 6-denominator topology (i. e. a graph in which all the propagators 
are different and all the numerators are equal to 1) related to the non-renormalized 
diagram of Fig. this topology is shown in Fig. El 

In the calculation, one also encounters the so-called subtopologies, corresponding 
to the topologies that can be obtained, from a given topology, by removing one 
or more propagator lines in all the possible ways. Starting from the topology of 
Fig. |2] and progressively removing one propagator, one finds the 4 independent 5- 
denominator topologies of Fig. El then the 6 independent 4-denominator topologies 
of Fig. EJ the 6 independent 3-denominator topologies of Fig. E] and finally the only 
2-denominator non- vanishing topology, the one shown in Fig. EJ 

In the graphical representation of subtopologies, internal straight lines corre- 
spond to propagators with mass m, internal wavy lines to massless (photon) propa- 
gators; external straight lines correspond to on mass-shell particles of mass m, while 
the momentum carried by external wavy lines is indicated explicitly. Note that the 
topology of Fig. Eland the topologies (c), (d) of Fig. El depend on both variables 
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Figure 2: The 6-denominator topology. 
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(a) 



(b) 



(d) 



Figure 3: The set of four 5-denominator topologies. The last one is the product of 
a 1-loop box and a tadpole. 



P 2 and Q 2 , the vertex subtopologies depend only on the square of the momentum 
written in the corresponding figures, while the subtopologies (e) of Fig. HI (d), (e), 
and (f) of Fig. and the one in Fig. [7|are constants (i.e. independent of P 2 , Q 2 ). 

The number of the subtopologies of any given topology can be large (as shown by 
the previous discussion); but different topologies can have common subtopologies, 
which amounts to say that many subtopologies can be known from independent 
previous work on graphs of different topologies. That is the case in the present 
calculation; as will be discussed later in more detail, most of the subtopologies of 
the current problem were indeed encountered and already worked out in JH] (which 
deals with vertex topologies). 



As it is well known ^T] the scalar integrals associated to any topology or 
subtopology are not all independent, as it is possible to establish several relations 
that link them among each other. As we will recall shortly in the following, there 
are (at least) three ways for writing such relations: using the Integration By Parts 
(IBPs) identities, exploiting the Lorentz structure of the integrals (LI) and relying 
on the symmetry properties (if any) of the integrals. 

We will use the following definition of the loop integration measure in D contin- 
uous dimensions, 



2.1 The Mis 
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Figure 6: The 2-denominator topology product of two 1-loop massive tadpoles. 



(corresponding to the energy scale /i = 1), where C(D) is the following function of 
the space-time dimension D: 

C{D) = {A*)^vU-^\ , (3) 

with the limiting value C(4) = 1 at D — 4. 

With this choice, the 1-loop tadpole with mass m reads 



k 2 + m 2 (D-2)(D-4) ' 1 1 

The most generic scalar integral associated to any 2-loop topology (or subtopol- 
ogy) can then be written as 
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(5) 



where r represents the number of different denominators D\, D T ; the z-th denom- 
inator is raised to the integer power rrii, with rrii > 1. In the numerator, there are 
q scalar products Si, S q , which involve one of the independent external momenta 
and one integration momentum kj, or two integration momenta. Since in the prob- 
lem at hand there are 3 independent external momenta and 2 integration momenta, 
there are 9 possible scalar products depending on the integration momenta; r of the 
scalar products can be simplified (or reduced) against the r propagators, so that the 
number of the "irreducible" scalar products remaining in the numerator is q = 9 — r. 

Having established the notation, we can describe the relations which hold for the 
integrals. 

• Integration by Parts Identities. 

Given any of the integrals defined in Eq. (jSj), its integrand can be used for 
writing the following 2 sets of identities [Hi] : 

f) ( s*™ 1 • • • S 1 ™ 9 1 

0%0% ^{_*_A_} = „, (6) 
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where the vector v M can be any of the 5 independent vectors of the problem: 
the momenta of integration k± and h^, and the external momenta pi, P2 and 
p%. In the context of dimensional regularization, every quantity appearing in 
Eqs. ()6I7|) is well defined, and the identities are always meaningful. Once the 
derivative acting on the integrands has been explicitly evaluated, Eqs. (|b17|) 
lead to a set of 10 identities for each initial integrand. These identities involve, 
as a rule, other integrals associated to the same topology as the initial one, 
where, at most, one of the powers of the scalar products and one of the 
powers rrij of the denominators can be one unit larger, while all the other 
powers remain the same or are decreased by one. In particular, it can happen 
that one of the propagators, raised to the first power in the original integrand, 
disappears as a consequence of the algebraic simplification against some re- 
ducible scalar product, generated by the differentiation; the resulting integrals 
are then associated to the subtopology where that propagator is missing. 

• Lorentz Invariance Identities. 

Another class of identities can be obtained from the fact that the integrals 
of Eq. (jSJ) are Lorentz scalars ^T] , and therefore they are invariant under in- 
finitesimal Lorentz transformations of the external momenta pi — > pi+5pi, with 
$Pi = e vPi> where the infinitesimal tensor is antisymmetric but otherwise 
arbitrary. That gives 



With the three independent external momenta of a four-point function, one 
can build three antisymmetric tensors of rank two; by saturating the above 
equation with the three tensors, one obtains the three identities 



Within the framework of dimensional regularization, it is possible in Eqs. (JO]- 
ITTj) to differentiate directly the integrand of the integral I{pi), before integra- 
tion. Eqs. (JUHTTj) then give three identities of a form similar to those obtained 
by the Integration by Parts method. 

• General Symmetry Relation Identities. 

Further identities among integrals can arise when the topology has some degree 
of symmetry. This can happen in cases in which, among the internal lines of 
a given topology, two or more represent particles of the same mass. In such 
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(j) (k) (1) (m) (n) 



Figure 7: The set of 14 master integrals involved in the calculation of the diagram 
of Fig. □ 



a case, there can be a transformation of the integration momenta which does 
not change the value of the integral, but changes the form of the integrand. By 
imposing the identity of the initial integral with the combination of integrals 
resulting from the transformation of the integration momenta, one can obtain 
additional identities relating integrals associated with the given topology and 
its subtopologies. An explicit example, involving the topology (b) in Fig. 0J is 
discussed in Section 2.0.3 of [TH] . 

For each topology, one can systematically write the above described identities 
starting from the integrand with all the powers rij of the scalar products equal to 
zero and all the powers of the denominators rrij equal to one, then in the case of all 
the integrands with N = Yli n i = 1 an d M = Ylj( m j — 1) = (with rrij > 0), then 
for all the integrands having N = and M = 1 (and always rrij > 0), then for all 
the integrands with N = 1 and M = 1 and so on. One finds that the number of the 
equations grows faster than the number of the involved integrals, until one obtains 
an apparently over-constrained set of linear equations for the integrals themselves. 
(For a more detailed discussion see for instance ^1]). The problem of solving such a 
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linear system (whose coefficients are polynomials in D, the masses and the external 
Mandelstam variables) is in principle trivial, but, due to the size of the system, 
it can be very demanding from the algebraic point of view (technical details on 
the computer programs developed in order to solve the system of linear equations 
are discussed in [IE])- As a result, one can identify a small number of so-called 
Master Integrals (Mis) for the considered problem, such that all the other integrals 
appearing in the considered identities are expressed as linear combinations of those 
Mis, with coefficients which are ratios of polynomials in D, masses and Mandelstam 
variables. It may also happen that all the integrals associated to a given topology 
can be expressed entirely in terms of the Mis of its subtopologies. 

For the purpose of this paper it is sufficient to say that all the scalar integrals 
that are necessary for the calculation of the 0(a 3 ) contributions to the Bhabha 
scattering amplitude from the considered Feynman graphs can be expressed in terms 
of 14 independent Mis only. There is some freedom in the choice of the integrals 
to be promoted to the role of Mis of the problem; we choose the set of Mis which 
are shown in Fig. [7| as " decorated graphs" . Each of the decorated graphs of Fig. [7| 
stands for a specific Master Integral; the denominator of the integrand is read from 
the lines (each line corresponding to a propagator raised to the first power, a line 
with a dot indicating that the corresponding propagator is squared), the numerator 
is given by the "decoration" ((p 3 • fc 2 ), (p3 • fci), etc. for graphs (b), (e), etc., or 
simply 1 when there are no other decorations). 

As a first remark, in Fig. [7|there is no 6-denominator MI: that means that all the 
scalar integrals associated to the 6-denominator topology of graph (c) in Fig. El can 
be expressed in terms of the Mis of its subtopologies with 5 or less denominators. 
Among the Mis of Fig. [3 the Mis (d)-(n) have already been calculated in 18 . As 
a consequence, in the present paper we focus our attention on graphs (a) and (b) 
of Fig. [7| Concerning graph (c) of Fig. it is the product of a massive tadpole, 
Eq. (j3J), and the 1-loop box graph. As the tadpole is singular as 1/(D— 4), one needs 
the 1-loop graph up to the first order in (D — 4) included; its complete calculation 
is provided in Appendix |Bj Let us observe that diagrams (a), (b), and (c) are the 
only Mis of the problem that depend on both the independent Mandelstam variables 
s = —P 2 and t = —Q 2 of Eq. (JTJ) (as well as on the electron squared mass m 2 ), while 
all the other graphs of Fig. [7| are functions of either P 2 or Q 2 (and m 2 ) only. 

3 The differential equations method 

The calculation of the Mis is performed by means of the differential equation method 
PHI |2~T) |23 El in this section, the main features of the method are recalled. 

Each of the Mis depends, in general, on all the independent Mandelstam variables 
of the problem, which we indicate by Sj = — p 2 , (i = 1,2,3,4), s 5 = — (pi — P3) 2 = 
—Q 2 , and s 6 = — (pi +P2) 2 = — P 2 (the invariants p 2 , p 2 , p 2 , and pi are constrained 
to be on the mass shell, p 2 = —m 2 , but that plays no role here). 
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When acting on any function of the Mandelstam variables, say M(s r ), one has 



Pi 



d 

M 



M(s r 



3 jp[ d Pk % 



:i2) 



where j,k = 1,2,3, while r = 1, • • • , 6. The ds^/dp^ are linear in the external 
momenta, so that the factors Pjds^/dp^ in the Eqs. (fT2"j) are linear in the Mandelstam 
variables. Eqs. (fT2|) can then be solved by expressing the d/ds^ in terms of the 

One obtains in this way 
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M{s r 



1 Q 2 + Am 2 

2 P 2 + Q 2 + Am 2 



Pi 
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dp^ 



P3 
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dps 



+ 



1 



P 2 



P 2 + Q 2 + Am 2 



P2 



d 



m 



Q 1 



d 



dQ- 



■M(s r ) 



P 2 + Q 2 +Am 2 

1 P 2 + Am 2 

2 P 2 + Q 2 + Am 2 

Q 2 



d 
dps 



Pi 



d 
dpj 



i- 1 - 

P2 



d 
dp^ 
d 



d 

dp2 



M(s r ), (13) 



+ 



P 2 + Q 2 + Am 2 



P2 



d 



m 



P 2 + Q 2 +Am 2 ^ 



Pi) 



d 
dp^ 



+ 



d 
dp^ 



+ 



d 

dp 11 



3/ J 



M(s r ). (14) 



We take one of the above equations, say Eq. (|T5jl for definiteness, and we replace 
the generic function M(s r ) by any of the Mis, say Mj(s r ). The l.h.s. is nothing but 
P 2 (dMi(s r )/dP 2 ); in the r.h.s., we write Mj(s r ) in its representation as an integral 
over the loop internal momenta, and we carry out the derivatives with respect to 
the momenta pj in the integrand. In this way, we obtain a combination of scalar 
integrals associated to the same topology as Mj(s r ), which, according to the previous 
discussion and the very definition of the Mis, can be expressed in terms of the Mis 
themselves. The result is a set of first order linear differential equations, of the form 



Q 2 -^M l (D,m 2 ,P 2 ,Q 2 



+ 



E 

3 

E 

k 



A% (D, m 2 , P 2 , Q 2 )M 3 (D, m 2 , P 2 , Q 2 ) 
(D, m 2 , P 2 , Q 2 )N k (D, m 2 , P 2 , Q 2 ) 

4? P, m 2 , P 2 , Q 2 )M 3 (D, m 2 , P 2 , Q 2 ) 

3 

E Bi S (D, m 2 , P 2 , Q 2 )N k (D, m 2 , P 2 , Q 2 ) 



(15) 



(16) 



where Mj(D,m 2 , P 2 ,Q 2 ) are the Mis of the topology under consideration, while 
N k (D,m 2 , P 2 ,Q 2 ) represent the Mis of the sub-topologies; finally, the coefficients 
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AfJ(D,m 2 ,P 2 ,Q 2 ) and B^(D,m 2 } P 2 ,Q 2 ) are ratios of polynomials in D,P 2 ,Q 2 
and m 2 . Note that the partial derivatives with respect to P 2 and Q 2 never mix 
within a same equation, so that the equations are in fact differential equations in 
a single variable. The two sets of differential equations in P 2 and Q 2 are therefore 
somehow redundant, and the redundancy can be used in the calculations as an a 
posteriori check. 

With some additional qualitative information on the Mis, the equations can also 
be exploited in order to fix the boundary conditions for the solution of the differential 
equations. We know that the considered box amplitudes are reg ular at P 2 = 0; 
therefore, by setting P 2 = in Eq. (fTtjj) (for arbitrary Q 2 ) the l.h.s. vanishes, while 
the r.h.s. gives a relation between the values of the considered Mj(s r ) at P 2 = and 
the values, at that same point, of the Mis of the subtopologies (which, in general, 
are simpler to obtain and supposedly known from previous calculations). 

While, in principle, the system of differential equations can be studied for arbi- 
trary values of the parameter D, we restrict our interest to the Laurent expansion 
of the Mis in powers of (D — 4). Therefore, we expand systematically in (D — 4) 
each of the Mis and of the coefficients that appear in Eqs. ()15I16|) . and solve the 
system directly for the Laurent coefficients of Mj(s r ). An explicit example of the 
procedure for the solution of the equations and the evaluation of the Mis is given in 
the following Section. 



3.1 The calculation of the 5-denominator box Mis 



The topology of Fig. |3] (c) is one of the topologies that present Mis which have not 
been already considered in ^H] (the second topology is the product of the 1-loop 
box, given in Appendix [0 and the tadpole). The two scalar integrals which have 
been chosen as the Mis for this topology are shown in Fig. [7| (a) and (b); according 
to the previous discussion, their explicit forms as loop integrals are 



F 2 (D,m 2 ,P 2 ,Q 2 




D D k{Z u k 2 



D, 



(pa -h)= / ^ D k^ D h 



V{DsD A VsD 6 

P3 ■ h 



V{D 3 V 4 V 5 V 6 



(17) 



By following the procedure outlined in the previous paragraphs, and dropping 
for ease of notation the arguments on which the two master integrals depend, one 
finds that the first order linear differential equations for Fi,F 2 in the variable P 2 
can be written as 

D-5 D-A 



op 2 

dF 2 

dP 2 



P 2 
1 



P 2 +Am 2 
D-5 



P 2 



+ 



P 2 + Q 2 +Am 2 
D-A 



Fx + ^A™ 2 ,^ 2 ,^ 2 ), (19) 
F 2 + n 2 (D,m 2 ,P 2 ,Q 2 ), (20) 



p 2 + 4m 2 p2 + Q2 +4m 2^ 

where fl\(D, m 2 , P 2 , Q 2 ) and fl 2 (D, m 2 , P 2 , Q 2 ) stand for lengthy combinations (not 
listed here for brevity) of the Mis corresponding to the graphs (c)-(n) of Fig. [7[ with 
coefficients given by ratios of polynomials in D,P 2 ,Q 2 m 2 . 
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The two equations of the system, Eqs. (|19I2(J|) . are completely decoupled, even 
for arbitrary value of D. The task of solving the system is then reduced to the 
integration of two independent first-order linear differential equations. The associ- 
ated homogeneous equation (which plays a key role in the solution of the equations) 
is exactly the same for both the Mis, a fact which further simplifies the explicit 
calculations. 

As already observed, the initial conditions are easily obtained by the equations 
themselves by imposing the analyticity of the solutions at P 2 = 0; in the case of 
Eqs. ()19I20|) . we can multiply both sides by P 2 and then take the P 2 = limit. The 
l.h.s. vanishes, while in the r.h.s., as the known terms Qi,Q,2 possess also a polar 
singularity 1/P 2 , we are left with (—1/2) times the values of the Fi at P 2 = and 
the residua of the singularities of the Qf, the explicit calculation gives 



F 1 (D,m 2 ,P 2 = 0,Q 2 ) = f\ \ , (21) 
F 2 (D,m 2 ,P 2 = 0,Q 2 ) = 




+ 



(5D-12) 1 (3D -7) 1 



m 2 2 (Q 2 + Am 2 ) 



3QD-2) 1 ST\ 

+ 2 m 2 (Q 2 + Am 2 ) K±J {kl ' ^ 



(3D -8) 1 



e 



32 m 2 

,[ 3(D-2) 2 1 1 1 Q m 

^64(D-3)m 4 1 , m 2 {Q 2 +Am 2 )) Q 

Let us observe that, as a check, one may obtain the two quantities Fi(D, m 2 , P 2 = 
0, Q 2 ) by setting pi = —p 2 (which implies P 2 = 0) directly in the definitions 
Eqs. fll9|2()j) of the two Mis. That leads to the following relations: 



F 1 (D,m 2 ,P 2 = 0,Q 2 )= f \ \ = I ^ D h^ D k 2 Pi p 3 p2p 6 . ( 23 ) 
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which is exactly Eq. (jUJ), and 



F 2 {D,m 2 ,P 2 = 0,Q 2 




V 1 V 3 V 2 5 V 6 



(24) 



where a dot on a propagator indicates that the corresponding denominator in the 
integrand is raised to the 2nd power. By expressing the above 4-denominator integral 
in terms of the Mis of Fig. [7[ Eq. (}2*2*|) is recovered. 

As pointed out in the previous section, we are interested in the Laurent expansion 
of the Mis F\ and F 2 with respect to (D — 4); it is known that they have at most 
double poles in (D — 4), so that 



Fi(P>, m 2 , P 2 , Q 2 ) = J]( J D-4) fc F 1 (fe) (m 2 ,P 2 ,Q 2 ) + 0( J D-4), (25) 

k=-2 

o 

F 2 (D,m 2 ,P 2 ,Q 2 ) = ^( J D-4) fc F 2 (fc) (m 2 ,P 2 ,Q 2 ) + G( J D-4). (26) 



k=-2 



By expanding in the same way also the inhomogeneous (known) terms, Eqs. (jl9l20j) 
generate a set of nested equations for the coefficients F- k \m 2 , P 2 , Q 2 ) of the expan- 
sion in (D — 4): 

9Ff } (m 2 ,P 2 ,Q 2 ) 



dP 2 



P 2 + P 2 + Am 2 



+^ k \m 2 ,P 2 ,Q 2 ). (27) 



where, due to P>-dependence of the homogeneous part of Eqs. (|19I2U|) . 



*l k \rn 2 ,P 2 ,Cf 



,(fc-i) 



(m 2 ,P 2 ,Q 2 



2 V.^ 2 + 4m 2 P 2 + Q 2 + Am 2 / 
+ fif } (m 2 ,P 2 ,Q 2 ) , (28) 

with Q^\m 2 , P 2 , Q 2 ) equal to the coefficient of order k in the Laurent-expansion of 
Qi(D, m 2 , P 2 ,Q 2 ) in powers of (D — 4); note that the complete inhomogeneous part 
of Eq. (|2Tj) contains (for k > —2) also terms in F[ k ^ l \m 2 , P 2 , Q 2 ) . 

The solution of the differential Eqs. (jl9|2())) . once written in the expanded form 
of Eq. (|2Tj) . is built, order by order in (D — 4), by repeatedly using Euler's method 
of the variation of the constants, see Eq. (|3T|) below. 

Euler's method requires the knowledge of the solution of the associated homo- 
geneous equation, which is the same for any order k of the expansion in (D — 4); in 
the case at hand it reads 

df(r) 



dr 



1 1 

r (r + 4m 2 ) 



f(r). 



(29) 
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whose solution, up to an irrelevant multiplicative constant, is 

f( r ) = 7=7 1 A ■>< • 

sj r[r + 4m z J 

The Euler's method then gives the solution of Eq. (J2*7)l in the form 



(30) 



F?\m 2 ,P 2 ,Q 2 ) 



v /p2(p 2 + 4m 2) 



/ dry/r{r + 4m 2 )^[ k) (m 2 , r, Q 2 ) 



(31) 



(k) 

the Kf J are the integration constants, which are fixed by imposing the initial con- 
ditions Eqs. pTI^j) at P 2 = 0. 

As we are interested in the expansion up to the finite part in (D — 4) and the 
expansion starts from 1/(D — 4) 2 , we need the first three terms of the expansion, 
i.e. we have to use repeatedly Eq. (JBTj) for k = —2, —1,0. It is actually convenient 
to replace the Mandelstam variables P 2 and Q 2 , by the dimensionless quantities x 
and y, defined as 



x 



VP 2 + Am 2 - VP* 



y 



VP 2 + 4:m 2 + VP 2 
and to introduce the functions 

F t {k \x,y) = Ft\m 2 ,P 2 ,Q 2 



y/Q 2 + 4m 2 - 
y/Q 2 + Am 2 + v/Q2 



(32) 



(33) 



The result can then be expressed in terms of 1- and 2-dimensional HPLs of argu- 
ment x and y and maximum weight 3; omitting for simplicity the dependence of the 
functions F^ k \x,y) on their arguments, and following the notation of Appendix IU1 
we have 



m 2 F ( - 2) 



m 2 F[~ x) 



1 

8 
1 

16 



1 



1 



1 - x) 
1 



(1 + x) 
1 



H(0;x), 



C(2) 



(34) 



m 2 F[ 0) 



_{l-x) (1+x) 
-H{0,0;x) + 2H(-l,0;x) 
1 1 



(1-1/) 



H(0;y) 



H(0;x) 



(35) 



1 

16 



C(2)+C(3)-21T(0;x)-C(2)H(-l;x) 



[l-x) (1+x) 
-H(0, 0; x) + 2#(-l, 0; x) - H(0, 0, 0; x) + H(-l, 0, 0; x) 
+H (0, 0; y) H(0; x) - 2#(-l, -1, 0; x) + H(0, -1, 0; x) 



(1-2/) 



4C(2)tf(0;y)+tf(0,0,0;y) 



+ 



C(2) -2tf(0; y) -4/7(0,0; y) -2/7(1,0; y)+6#(-l,0; y))#(0; ar) 
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,(-2) 



<0) 



+ (3C(2) +#(0,0; y)j (G(-y; x)-G(-l/y; x)j + 2H(0; y)H(-l,0; x) 
-H(0; y) (G(-y, 0; x) + G(-l/y, 0; x)) +G(-y, 0, 0; x) 

-G(-l/y,0,0;x)]\, (36) 
' ' ' (37) 



1 

32 
1 

64 



(1-x) (1 + x 
1 1 



y 



2 + y 



_(l-x) (l+x)_ 
+2H(-l,0;x 



1 

~-y 

y 



H(0;x), 

C(2)-4tf(0;x)-tf(0,0;x) 
H(0;y)H(0;x) 



--2+y 

y 



+ 



32 



i <2 



2 - 



(1-x) 
1 



(1-x) 



H(0;x), 
C(2) - 2 



(i-y) 0- + v) 



+H(0, 0, 0; y) + H(0, 0; y) - 5H(0; x) + 2H(-1, 0; x) 



(1-x) 
1 

(1 +x) 
1 



5H(0; x) + H(0, 0; x) - 2H(-1, 0; x) 



+ 



(1 - x) 



(C(2) - 2F(0, 0; y))#(0; x) + #(0, 0, 0; x) 



128 



(1-x) (1 + x) 



--2 + y 
.2/ 



4C(2) + 2C(3) 



-2C(2)F(-1; x) - 2(7 - H(0, 0; y))ff(0; x) - 4iJ(0, 0; x) 
+8H(-1, 0; x)-2#(0, 0, 0; x) + 2H(-l, 0, 0; x)+2H(0 -1, 0; x) 





"1 




+ 


- -y 
.y 





4((2)H(0;y) + H(0,0,0;y) 



-4#(-l,-l,0;x) 



-2(2/7(0; y) - -C(2) + 2#(0, 0; </) + H(l, 0; </) 

-3if (-1, 0; y))H(0; x) + (3((2) + H(0, 0; x) 
-G(-l/y; x)) + 2#(0; y)if (-1, 0; x) - H(0; y)G(-y, 0; x) 

-H(0; y)G(-l/y, 0; x)+G(-y, 0, 0; x )-G(-l/y, 0, 0; x) 



(38) 



(39) 



As a check of the results reported above, the expressions for the Mis F[ k) (x,y) 
and F^\x,y) have been inserted in the corresponding differential equations with 
respect to the variable Q 2 (Eq. HK|) : those equations were found to be satisfied. 

The results for the Mis can be downloaded as an input file for FORM in j33j . 



15 



3.2 Asymptotic expansions 



In this Subsection we provide the asymptotic expansions, in various kinematic re- 



gions, of the coefficients and F 2 (K> ( k = —2, —1,0) introduced above. 

The first region of interest is the one in which P 2 ~ Q 2 ^> m 2 , (i ~ j/ < 1), 
relevant for the large angle Bhabha scattering. Employing the definitions 

and keeping only the leading terms, we find 



, / m 2 
In 



m' 
Q 2 " 



m 2 F[~ 



m 2 F (o) 

1 96 P 2 



,(-2) 



,(0) 



1 m 2 

1 — [2C(2) 
16 P 2 L 



4L r — L 2 — 2L r L w \ j , 



' ;;r ' r i2C(2) + 12C(3) + 18C(2)lnfl + ^) + G«2)Z, 



-24L r - 6L 2 + 3L In 1 + 



-C)L r L, r hi ( 1 - ^ 



Q 



P 2 

6L r L 2 +6L r Li 2 



P 2 



— 2L — \2L r L w 



-3/4 In (J-^ ) / - GI (r I/ 2 

16 P 2 r ' 

32 r + 64 P 2 
1 

128 



P 2 ^ 



P 2 



P 2 

+ 6Li 3 



4C(2) + 10L r + 2L 2 + L 



2£(2) — 8L r — L 2 — 2L r L w 
l Q 2 ~ 



128 P 2 



8C(2) + 4C(3) 

Q 2 



+6C(2)ln( 1 + ^ )-14L r + 2C(2)L r -4L 2 +L 2 ln ( 1 + 



P 2 



2 L 3 

3 r 



lL r L w - 



■2L r L w \n I 1 + ^7 



+8C(2)L U , + L 2 In 1 + 



P 2 



Q 2 
P 2 

l 



■2L r L 2 w + 2L r Li 2 

Q 2 



P 2 



+ -L„, - 2L w Li 2 - 



P 2 



+2Li 3 - 



P 2 



(40) 



(41) 
(42) 



(43) 
(44) 
(45) 



(46) 



The case P 2 3> Q 2 3> m 2 can be immediately obtained from the previous equa- 
tions. 

The second region of interest is the one in which again P 2 ^> m 2 , while Q 2 is 
much smaller than P 2 (Q 2 <C P 2 ), but otherwise arbitrary, so that i « i/ < I. 
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Keeping only the leading terms we find: 



m 2 Ft 2) 



m 2 F[- 1] 



m 2 F± 0) 



,(-2) 



,(0) 



1 m 2 
4^ n 



1 m 2 
16 P 2 

1 m 2 
~16P2 



2C(2) - 2 



4m 2 



tf(0;y) 



(47) 
(48) 



2C(2) + 2C(3) + ^/l + ^- 2r 



4C(2)if(0; y)+if (0,0,0;y) 



+ 



/ Ain 2 / 

4 - 2H(Q, Q ]y )-^l + — [C(2) - 2H(0; y) - AH(0, 0; y) 

1 



+6if(-l,0;j/)-2if(l,0;y)) 



16 P 2 r 



(49) 
(50) 



Q1 f l + -^)H(0;y)LA, 



Am 2 



Am 2 



(51) 



1 

64 



2((2) + Pf (0,0;y) + 



+H(o,o,Q;y) +5L r + L 2 
1 



. 1 + 

1 m 2 
32 P 2 



-2H(0,0;v) 



Q 2 



L r — -zL — -L H - 

2 6 m 2 



4 + C(2)- [2 + C(2) 
2C(2) + C(3) - (7 



-tf(0,0;y))L r -L 2 --L r 3 



^( 1 + 4S)I 8 »^ 
+2Pf (0, 0, 0; y) + (2((2) -8if (0; y) - 8tf (0, 0; y) 

-4//(l,0;y) + 12i/(-l,0;y))L r 
In the extreme case y — > 1, i.e. m 2 > Q 2 > 0, the previous expansion gives: 



m 2 F^ 


1 m 2 

= IT 2 "^ 


m 2 ^ 


1 m 2 " 
16 P 2 . 


m 2 F (1) 


1 m 2 ■ 
48 P 2 . 




= o, 



2C(2)-L 2 j 

18C(2) - 6C(3) + 12L r + 3F 2 + L 



(52) 

(53) 

(54) 

(55) 
(56) 
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F (2) 



,(2) 



L r + 

32 



1 m 
16 P 2 



1 - L, 



1 

64 L 



4C(2) + 5L r + L 2 
l Q 



-3Li - L 



+ 



384 m 2 



1 m 2 
" 192 P 2 

5-2C(2) 



24 + 6C(2) - 12L r -6C(2)L r 



(57) 



(5* 



4 The scalar 6-denominator integral 

We have already observed that there is no MI associated to the 6-denominator 
topology (c) in Fig. |21 The scalar integral associated to the original graph in Fig. ^ 
(which has 6 propagators, one of them raised to the second power), therefore, is not 
a MI and can be expressed in terms of the 5, 4, 3 and 2-denominator Mis of Fig. 
For completeness, we report its definition and analytic value: 



1 



J2(D-4) k AW(x,y)+0((D-4) 



(59) 
(60) 



k=-2 



where the coefficients A^ k \x,y) of the Laurent expansion in powers of (D — 4) are 
given (dropping again for simplicity the arguments x, y) by: 



m 



m 



6^(-2) 



4(1-?/) 2 



8(l + z)(l-y) 
1 



(1-x) (1 + x) 
1 



(l-yf 

1 



H(0;x) 



(61) 



(1 + x) 



1 - 



(i-y) 



16(l + x) 2 (l -y) 



(l + x)_ 



1 



(1-2/) 



8(1 -y)3 
1 

48(l-x)(l-y) 
1 



_(l-x) (1+;X 
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2 - 



1 



16(l-2/) 2 
1 



(1-2/) 
1 



(1-2/) 
1 



1 + 



(l+x)_ 
1 1 



m 



> A (°) 



4(1 -y) 2 
1 

(l 3 ^ 
1 



(1-x) (l + x) 

fl(i;i/)W;i), 



24(l+ar)(l-y) 

3 2 
(l_ y ) 2 + (l_ y )3 

2 



(1-2/) 
H(0;y) 



(1-2/) 



(62) 



(l-y)L(l-x) (x+l)J|_(l-y) (l-y) 2 \{ 48 



r #(0;x)if(0;j/) 
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(1 _ y)2 j L ff(0;y)(2C(2)+ff(-l,0;x)-^(0;x) 
-^if (0; y) (G(-l/y, 0; x) + G(-y, 0; x) + ^H(0; x) + 4((2) 



+ 



1 

+ 8 



1 



(1-x) (x+1) 

1 1 

(l-y)~(l-yy 



(1-y) (l-y) ; 



-^(3C(2)+ff(0,0;x; 



1 



G(-l/y, 0; x) +G(-y, 0; x) + -if (0; x) 



+4C(2))( i //(0;y) + //(l; 2/ ) 



1 

48 



1 

16 



x 



+ 



(i-y) (i-y) 2 
l 



l 



i + 



x + 1)_ 

1 1 



(x + 1) (x+1) 2 
2 



(x + 1) 2 



+ 6 



3C(2) + #(0,0;x)^ 

3C(2)+tf(0,0;x) 
1 



x 



.(1-x) (x+1) 





1 1 






[(1-y) (1-y) 2 . 


|24 



(1-y) (1-y) 2 . 

//(-l,0;x)+2C(2) 



-^ff(0;x) 



1 

16 



1 



.(1-y) (1-y) 2 
+tf(0;y)+4tf(l;y))4 
16 



- ^^(0; x) (l + 4/f (0; y) + 4/7(1; y)) 

1 1 r ^2C(2) + /J(-l,0;x)-^(0;x)^(2 

ff(0;x)(-6ff(-l,0;y)-6C(2) 



+ 



16(l-y) 

yif (0; y) + 4/7(0, 0; y) + 2ff (1, 0; y)) - 4C(2)ff(0; y) 

-if (0, 0; y) (G(-y; x) - G(-l/y; x)) - if (0, 0, 0; y) 
-3C(2)G(-y; x) - G(-y, 0, 0; x) + 3C(2)G(-l/y; x) 

+G(-l/y, 0, 0; x)j + ^ (^(3) - 6C(2)ff (-1; x) 

3 



-12ff(-l, -1, 0; x) + if (0; x)\ll- -C(2) + 9/f (-1, 0; y) 
-8ff (0; y) + 3ff (1, 0; y) + 12ff (1, 1; y) + 6ff (0, 1; y)ff (0; x) 
-3ff(0, 0; y)ff (0; x)) + 6ff (-1, 0, 0; x) - 2C(2)ff(0; y) 
3 



+ 



-G(-y; x) (3C(2) + if (0, 0; y) + 4/7(0, 1; y)) 

G(-l/y; x) (7C(2) +ff (0, 0; y) +4/f (0, 1; y)) - |if (0, 0, 0; y) 
-6if (0, 0, 1; y) + 6G(-y, -1, 0; x) - ^G(-y, 0, 0; x) 



+6G(-l/y, -1, 0; x) - ^G(-l/y, 0, 0; x)) | 
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1 



(1-y) (l-y) 2 \ (z+l) 



1 

24 



1 

16 



(x + 1) 



2({2)--H(0;x)+H(-l,0;x) 
6 



(63) 



5 Summary 



The aim of the present paper was to identify and evaluate the Mis which are nec- 
essary for the calculation of the contribution of two-loop box diagrams with one 
electron loop (Np = 1) to the corrections of 0(a 3 ) to the Bhabha scattering am- 
plitude in QED. The result was obtained without any approximation, so that our 
expressions keep the full dependence of the Mis on the electron mass m and on the 
Mandelstam variables s,t. 

It has been shown that all the scalar integrals occurring in the problem can be 
expressed in terms of a set of 14 Mis; the reduction procedure has been carried out 
with the by now standard approach based on the use of the IBPs, Lis, and symmetry 
identities. 

Out of the 14 Mis of the set, 11 had already been calculated in a previous work. 
For what concerns the remaining three, one is simply the product of two 1-loop 
integrals, while the other two are genuine 2-loop integrals; they all depend on both 
the Mandelstam variables s and t. 

The central part of the paper has been devoted to the evaluation of the two Mis 
not already known in the literature. The calculation has been performed with the 
method of the differential equations with respect to the Mandelstam variables. While 
a two by two system of differential equations was expected on general grounds for 
the two Mis, it was found that the two equations of the system are in fact decoupled, 
greatly simplifying the task of finding the solution. The boundary conditions for 
the solutions of the differential equations corresponding to the considered Feynman 
graph integrals have been explicitly obtained from the equations themselves and the 
qualitative knowledge of the analytical properties of the Mis, namely the regularity 
of the Mis at P 2 = 0. 

The analytic expression for the Mis has been given as a Laurent series in powers 
of (D — 4), where D is the dimensional- regulator for both IR and UV divergences. 
The coefficients of the Laurent series of order \/{D — 4) 2 , \/{D — 4), and zeroth order 
in (D— 4), have been written in closed analytic form in terms of 1- and 2-dimensional 
harmonic polylogarithms, of maximum weight w — 3. 

For completeness, the Laurent expansion for the scalar 6-denominator integral 
associated to the graph in Fig. ^ is also given up to the same order in [D — 4). 

With these results, it is now possible to evaluate the contributions of the two-loop 
box diagrams with one electron loop to the amplitude of the Bhabha scattering in 
QED, at 0(a 3 ). The explicit expression for the 0(a 3 ) contribution to the amplitude 
will be given elsewhere. 
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A Propagators 

In this Appendix we list the dependence on the loop momenta of the denominators 
of the propagators appearing in the various integrals. 



T>x = kj, (64) 

V 2 = (Pi-ps-h) 2 , (65) 

V 3 = [k 2 2 +m 2 ], (66) 

V A = [(p! - h) 2 + m 2 ] , (67) 

A, = [(P2 + h) 2 + m 2 ] , (68) 

= [(pi - P3 ~ h + k 2 f + m 2 } . (69) 



B One-loop results 

Among the 14 Mis of Fig. [7J there are 4 Mis which are simply products of two one- 
loop diagrams: the MI shown in Fig. [7| (c) is the product of a box diagram and a 
tadpole; the graph (h) is the product of a vertex in the t-channel and a tadpole; the 
one in (k) is the product of a tadpole and a bubble in the t-channel with massless 
internal lines. Finally (1) is the product of a tadpole and a a two-point function in 
the s-channel with massive propagators. 

While the value of the tadpole was given in Eq. (jij), we refer to ^H] for the 
explicit expressions of the massless and massive bubble diagrams and the vertex 
correction. We discuss here the calculation of the one-loop box diagram, up to the 
first order in (D— 4) included. The considered topology (which in this case coincides 
with the MI itself) is shown in Fig. |HJ 

By applying the reduction procedure outlined in Section it is found (as ex- 
pected) that the only MI for that topology is the scalar integral B(D, m 2 , P 2 , Q 2 ): 
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Figure 8: The 1-loop box. 



where, as usual, P 2 = (pi + p 2 ) 2 and Q 2 = (pi — P3) 2 , such that P 2 = —s and 
Q 2 = -t. 

By following the method outlined in Sectional we find that B(D,m 2 , P 2 ,Q 2 ) 
obeys the following first-order linear differential equations in P 2 and Q 2 : 



dB 

W 2 

dB 
W 



(D-5) 



(D-4) 



P 2 (P 2 + 4m 2 ) (P 2 + Q 2 +4m 2 ) 
(D - 6) (D — 4) 



Q 2 



(P 2 +Q 2 +4m 2 ) 



B + n^D.m 2 ^ 2 ^ 2 ), (71) 
B + n 2 (D,m 2 ,P 2 ,Q 2 ), (72) 



where the homogeneous part of the differential equation in P 2 is the same as in 
Eq. and the non-homogeneous terms fli and Q2 (not to be confused with the 

inhomogeneous terms of Eqs. ()19l20j) !) are given by 



£li(D, m 2 , P 2 , Q 2 ^ 



-{D-A) 



1 



(Q 2 + 4m 2 



_4m 2 P 2 Am 2 Q 2 (P 2 + Am 2 ) 
1 



+ 



Q 2 (P 2 + Q 2 + 4m 2 ) 




2(P>-3) 



Q 1 



1 



Q 2 (P 2 + Q 2 + 4m 2 ) 



(P 2 + 4m 2 ) 2 Q 2 (P 2 + 4m 2 ; 

O 



+ 



P-3) 
2m 2 Q 2 

m 2 Q 2 



1 



P 2 (P 2 + 4m 2 ) 
1 




(P 2 + 4m 2 ) 2 Q 2 (P 2 + 4m 2 ) 



1 



Q 2 (P 2 + Q 2 + 4m 2 ) 



(73) 
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n 2 (D jm \p 2 ,Q 2 " 



(D-A) 



(P 2 + Am 2 ) [Q 2 (P 2 + Q 2 + Am 2 )_ 



2(D-3) 




(P 2 + Am 2 ) 2 



_Q 2 (P 2 + Q 2 + Am 2 )_ 



o 



(D-2) 



m 2 (P 2 +Am 2 ) 2 [Q 2 (P 2 + Q 2 +Am 2 )_ 



(74) 



Each of the two equations Eqs. ()71|72j) is sufficient to obtain the explicit value 
of B(D, m 2 , P 2 , Q 2 )] we use the equation in P 2 for definiteness. Following the lines 
of Section Kill we find as boundary condition at P 2 = 



B(D,m 2 ,P 2 = 0,Q 5 



2m 2 




(D-3) 



m 2 Q 2 




(75) 



The differential equation Eq. (J71j) can be explicitly solved by the method of 
the variation of the constants of Euler outlined in Section 13.11 The homogeneous 
equation, as already observed, is the same as in Eq. (fT^j). so that the coefficients of 
the Laurent expansion in (D — A) of B(D, m 2 , P 2 ,Q 2 ) can be obtained as in Eq. (J3T| . 

In terms of the dimensionless variables x and y defined in Eq. ()32j) . the Laurent 
expansion in (D — A), which begins with a single pole, is 



J2(D-A) k B^(x,y) + 0((D-A) 2 ) . 



(76) 



Since the one-loop scalar box graph is multiplied by a tadpole, which is singular as 
l/(D — 4), in the MI of Fig. (c), we need its Laurent expansion up to the first 
order in (D — A) included. 

Dropping for simplicity the dependence on (x, y) of the coefficients of the Laurent 
expansion B^(x,y) we find the following results 



1 



2 L(l-2/) 



1 



1 



L(l-y) (1-y) 2 
H(0;y))H(0;x), 



[l-x) 
1 



(l + x) 

i 

JTxj 



m 



4 S (1) 



H(0;x), 
-2H(l;y) 

-2C(3) 



(77) 



(78) 



(1-y) (l-y) 2 \ [(l-x) (l+x) J 
+2C(2)#(-1; x) + AH(-1, -1, 0; x) 
+2H(-l,0;x)[H(0;y) + 2H(l;y)} 

-2H(-1, 0, 0; x) + H(0; x) [((2) - 2H(1, 0; y) - AH(1, 1; y)\ 
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+H(0; y) [4C(2) - G(-y, 0; x) - G(-l/y, 0; x)) 

- [H(0, 0; y) + 2H(0, 1; y)] [if(0; x) - G(-y; x) + G(-l/y; x)} 

+H(0,0,0;y)+2H{0,0,l;y) 

-2H(1- y) [G{-y, 0; x) + 0; x)) + 3((2)G(- y] x) 

-2G{-y, -1, 0; x) + G{-y, 0, 0; x) - 5C(2)G(-l/y; x) 

-2G(-l/y, -1, 0; a;) + G(-l/y, 0, 0; x) } . (79) 

Eqs. (l77fJT9~j) are valid in the non-physical region P 2 = — s > 0; the corresponding 
expressions for the physical region are recovered by standard analytical continuation. 



C Harmonic Poly logarithms 

In this Appendix we briefly review some of the properties of the Harmonic Polylog- 
arithms of one variable, x, (HPLs), introduced in [21] as an extension of Nielsen's 
polylogarithms (2U HUH OHl IH2] , as well as of the Harmonic Polylogarithms of two 
variables x and y (2dHPLs), introduced in fTEL 



C.l One-dimensional Harmonic Polylogarithms 

One starts by defining the following set of algebraic factors 

/( - 1:x) = (TTi) ' (80) 

f(0-,x) = -, (81) 

X 

/&;«) = j^-y (82) 

The one-dimensional HPL, H(m w ; x), can then be defined as the set of functions 
generated by the repeated integrations 



f dz {/(-l; z); /(0; /(l; *)} # K„; z) 
Jo 



(83) 



with 



= / = ln(l + x), (84) 



o 

H(0;x) = lux, (85) 

H(l;x) = [ X j-^- = -ln(l-s), (86) 
Jo 

and where is a w-dimensional vector whose components can assume the values 
1, or —1; w is said the weight of the corresponding HPL. 



24 



The following relations are valid: 



H(O w ,x) = -^ln w x, (87) 

ex 



\ dzf(a;z)H(m w _ 1 ;z) 
Jo 



^-H(a ) m w _ 1 ;x) = f(a;x)H(m w _ 1 ;x) , (89) 
ax 

where W stands for the vector whose w components are all equal to 0. 

The set of the HPLs fulfills an algebra; the product of two HPLs of the same 
argument x and of weights W\ and W2 is a suitable combination of HPLs of the same 
argument and weight w = w\ + W2- 

H(p;x)H(q;x)= ^ H(r; x) , (90) 

r=p+q 

where r is a (w p + w g )-dimensional vector constituted by all mergers of p and q 
in which the relative orders of the elements of p and q are preserved. In the case 
w p = 1, for example, we have: 

H(a;x)H(mi, ■ ■ ■ ,m q ;x) = H(a, mi, ■ ■ ■ ,m q ; x) + H(mi, a, rri2, ■ ■ ■ ,m q ;x) 

H \-H(mi, ■ • • ,rn q , a;x) . (91) 

Moreover, by subsequent integrations by parts on the definition itself, the fol- 
lowing identities between HPLs hold: 

H(mi,---,m q ;x) = H(mi,x)H(m2, ■ ■ ■ ,m q ;x) 

-H(m 2 , mi, x)H(m 3 , ■ ■ • , m q ; x) 

-\-H(m3, m2, mi; x)H(m4, ■ ■ ■ 

(-l)^(m 9; ..., mi ;x), (92) 

For a more complete treatment (in particular for the analytical continuation) 
and the numerical evaluation of the HPLs we refer the reader to EH and l25l. 



C.2 Two-dimensional Harmonic Poly logarithms 

We recall here shortly the definition and properties of the 2 dimensional HPLs 
already used in Eqs. (j36l39j) . As observed in they can be simply obtained by 
replacing the f(i; x) of Eqs. (jHOHHH) by a generalized set of factors 

g(i;x) 



x — % 



where in the present calculation the "index" i spans the enlarged set of values 0, 
-1, -y and -1/y: 

»(-!;*) = J^y 03) 
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9(0; x) = -, (94) 

x 

g(-y;x) = — L_, (95) 

9(-l/v,x) = -7——rr- (96) 

2dHPLs are then denned as the set of functions generated by the repeated inte- 
grations 

px 

(97) 



px 

I dz{g(j;z)} G(m w ;z) 
Jo 



where j and the components of m„, can take the values 0, —1, —y, and —1/y. In 
particular 

G(-l;x) = \n(l + x) = H(-l;x), (98) 
£7(0; ar) = \nx = H(0;x), (99) 

G(-y;x) = [ X = \n(l + -), (100) 

Jo (z + y) V y) 

G(-l/y;x) = / s; = ]n(l + xy). (101) 

It is not difficult to see that the 2dHPLs involving the subset of indices 0, (—1/y) , 
(— 1) and argument x can be re-expressed in terms of ldHPLs and 2dHPLs of argu- 
ment xy and indices 0, (—y) and (—1); in particular we have 

G(-l/y;x) = H(-l;xy), (102) 

G(-l/y,0;x) = H{-1, 0;xy) - H{0;y)H{-l; xy) , (103) 
G(-l/y,0,0;x) = H{-1, 0, 0; xy) + H(0, 0;y)H(-l; xy) 

-H{0;y)H{-l,0;xy), (104) 

G{-l/y,-l,0;x) = G{-1, -y, 0; xy) - H{0;y)G{-l, -y; xy) . (105) 

The analytical continuation of the 2dHPLs listed above can be obtained by following 
the lines of [21! • 

Concerning their numerical evaluation, the results of |2Sj cannot be used here 
as they apply to a different set of indices of the 2dHPLs. For that reason, we give 
now their expressions in terms of Nielsen's polylogarithms of non-trivial argument. 
Let us remind here that such expressions are by no means unique, as those polylog- 
arithms of non-trivial argument can satisfy several identities, often quite involved; 
an elementary example is 

Li 2 (l - z) + Li 2 (z) + log z log (1 - z) - C(2) = 0. 

The expressions in terms of 2dHPLs of suitable indices and argument x do not suffer 
of that drawback. 
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We find 



G(-y,0;x) 
G(-l/y,0;x) 
G(-y, 0,0; x) 



G(-l/y,0,0;x) 



G(-y, -1,0; x) 



ln(s)ln(l + ^+Li 2 (- X - 
In (x) In (1 + xy) + Li 2 (— xy), 



dx' 1 



InV 



x' + y 2 



^,( 1 + £) +lnWLi2 (_£)_ Li3 (_£ 



<ix' 1 



InV 



o + 1 2 
ln 2 x 



2 

* dx' 



In (1 + xy) + lnxLi 2 (— xy) — Li 3 (— xy), 

x ' dr" 

—r. r lnx 

o % + y Jo x" + i 

ln(l+x) In (: 



(106) 
(107) 
(108) 

(109) 

(110) 

(111) 
(112) 



ln(x) + ^ln(l+x)-^ln + -Info) 
+ ln(l + x) Info) - \ ln(x) ln(l + x) - \ ln 2 (l + x) 
-i ln(l+x) ln(l-y) + iln 2 (l-y)-ln(l-y) Info) 



x 

V ' y. 



i 



+ 



ln(l+x) 

\ ln(x) In + +i In (l + ^) ln(l - y)-i ln 2 (l - y) 
^ In ^1 + ^ Info) +ln(l-y) Info) In (\ + 



ln(x) - In ( 1 + - ) + ln(l - y) - Info) 



Lio(-x) 



+ 



ln(x)+ln(l-y)-lnfo) 



Li 2 (--)-ln(l + -)Li 2i 

yj \ yj \%+y. 



y 



-lnfo)Li 2 fo) - 



ln(l + x) - In 1 + 



ln(l + x) - In 1 + 



x 



x 



Info) 



Li 2 



Li 2 



x + y 



x + y 



1 + x 



+ 



ln(x) — ln(l + x) + ln(l — y) — Info) 



Li 2 



s(i -y) 

x + y 



+Li3 (- I )- Li 3(-£) + Li 3 (-^) + Li3 W 



-Li 3 

+Sl, 2 



x + y 
x + y 
1 + x 



+ Li 3 



(1 +x)y 
x + y 



x + y 
1 + x 



Si, 2 (y) 



(113) 
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G(-l/y-l,0;x) 



dx' 



dx" 



x' + 1 Jn x" + 1 



lnx" 



(114) 



= -C(3) + 



Ml/) 



ln(x) ln(l + x) - - ln 2 (l + x) + - ln(l + x) ln(y) 



ln(l - y) \n(y) 
1 



ln(l + x) + 



2 ln(x) ln(l + x) - ln 2 (l + x) 



ln(x) ln(l + xy) + ln(l + x) ln(l + xy) 



--ln 2 (l + xy) 



ln(l + xy) + ln(x) +ln(l -y)]Li 2 (-xy) 
ln(x) + ln(l — y) — ln(l + xy) 
ln(x) — ln(l + x) + ln(l — y) 



Li 2 f-x) 



Li 2 



x(l - y) 



1 + x 



ln(l + x) - ln(l - y) 



Li 2 



(1 + x)y 



+ 



ln(l +x)-ln(l + xy) 
x(l - y) 



+Li 3 (-x) + Li 3 



1 + x 



Li 3 (-xy) 



Li f 



1 + xy 

1-2/ 
(1 + x)y 



-Li a 



(l+x)y' 
1+xy 



') +Li3 (tS) + S i 5 2(2/) + Si, 2 (-^). (115) 



D Limiting cases 



In order to be able to impose the initial conditions for the solutions of the differential 
equations for the Mis, we need the expressions of the 2dHPLs in some particular 
point of one of the variables. In our case we needed the following values of 2dHPLs 
at x = 1 in terms of HPLs of argument y: 



G(-y; 1 
G(-l/y;l 
G(-y,0;l 
G(-l/y,0;l 
G(-y,0,0;l 
G(-l/y,0,0;l 

G(-y,-l,0;l 



H(-l;y)-H(0;y), 
H(-l;y), 

-C(2) + H(0,-l;y)-H(0,0;y), 
-#(0,-1; y), 



(116) 
(117) 
(118) 
(119) 



C(2)tf(0; y) - H(0, 0, -1; y) + H(0, 0, 0; y), (120) 

#(0,0, -l;y), (121) 

-^(3)+ 1 -C(2)[H(l;y)-H(-l;y)\ 

-#(i,o,-i;y) + #(i,o,0;y), (122) 
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£(-!/</, -1,0;1) = ~C(2)[fT(l;y) + J ff(-l;y)] + J ff(0,0,-l;y) 



+#(1,0,-1;!,). (123) 
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